Predicting the Long-Term Effects of Therapeutic Neutralization of Oncostatin M on Human Hematopoiesis.

Therapeutic neutralization of Oncostatin M (OSM) causes mechanism-driven anemia and thrombocytopenia, which narrows the therapeutic window complicating the selection of doses (and dosing intervals) that optimize efficacy and safety. We utilized clinical data from studies of an anti-OSM monoclonal antibody (GSK2330811) in healthy volunteers (n = 49) and systemic sclerosis patients (n = 35), to quantitatively determine the link between OSM and alterations in red blood cell (RBC) and platelet production. Longitudinal changes in hematopoietic variables (including RBCs, reticulocytes, platelets, erythropoietin, and thrombopoietin) were linked in a physiology-based model, to capture the long-term effects and variability of therapeutic OSM neutralization on human hematopoiesis. Free serum OSM stimulated precursor cell production through sigmoidal relations, with higher maximum suppression (Imax ) and OSM concentration for 50% suppression (IC50 ) for platelets (89.1% [95% confidence interval: 83.4-93.0], 6.03 pg/mL [4.41-8.26]) than RBCs (57.0% [49.7-64.0], 2.93 pg/mL [2.55-3.36]). Reduction in hemoglobin and platelets increased erythro- and thrombopoietin, respectively, prompting reticulocytosis and (partially) alleviating OSM-restricted hematopoiesis. The physiology-based model was substantiated by preclinical data and utilized in exploration of once-weekly or every other week dosing regimens. Predictions revealed an (for the indication) unacceptable occurrence of grade 2 (67% [58-76], 29% [20-38]) and grade 3 (17% [10-25], 3% [0-7]) anemias, with limited thrombocytopenia. Individual extent of RBC precursor modulation was moderately correlated to skin mRNA gene expression changes. The physiological basis and consideration of interplay among hematopoietic variables makes the model generalizable to other drug and nondrug scenarios, with adaptations for patient populations, diseases, and therapeutics that modulate hematopoiesis or exhibit risk of anemia and/or thrombocytopenia.

Quantitative modeling of tumor dynamics and development of drug resistance in non-small cell lung cancer patients treated with erlotinib.

Insight into the development of treatment resistance can support the optimization of anticancer treatments. This study aims to characterize the tumor dynamics and development of drug resistance in patients with non-small cell lung cancer treated with erlotinib, and investigate the relationship between baseline circulating tumor DNA (ctDNA) data and tumor dynamics. Data obtained for the analysis included (1) intensively sampled erlotinib concentrations from 29 patients from two previous pharmacokinetic (PK) studies, and (2) tumor sizes, ctDNA measurements, and sparsely sampled erlotinib concentrations from 18 patients from the START-TKI study. A two-compartment population PK model was first developed which well-described the PK data. The PK model was subsequently applied to investigate the exposure-tumor dynamics relationship. To characterize the tumor dynamics, models accounting for intra-tumor heterogeneity and acquired resistance with or without primary resistance were investigated. Eventually, the model assumed acquired resistance only resulted in an adequate fit. Additionally, models with or without exposure-dependent treatment effect were explored, and no significant exposure-response relationship for erlotinib was identified within the observed exposure range. Subsequently, the correlation of baseline ctDNA data on EGFR and TP53 variants with tumor dynamics' parameters was explored. The analysis indicated that higher baseline plasma EGFR mutation levels correlated with increased tumor growth rates, and the inclusion of ctDNA measurements improved model fit. This result suggests that quantitative ctDNA measurements at baseline have the potential to be a predictor of anticancer treatment response. The developed model can potentially be applied to design optimal treatment regimens that better overcome resistance.

Increase in peg-asparaginase clearance as a predictor for inactivation in patients with acute lymphoblastic leukemia.

Asparaginase is an essential component of acute lymphoblastic leukemia (ALL) therapy, yet its associated toxicities often lead to treatment discontinuation, increasing the risk of relapse. Hypersensitivity reactions include clinical allergies, silent inactivation, or allergy-like responses. We hypothesized that even moderate increases in asparaginase clearance are related to later inactivation. We therefore explored mandatory monitoring of asparaginase enzyme activity (AEA) in patients with ALL aged 1-45 years treated according to the ALLTogether pilot protocol in the Nordic and Baltic countries to relate mean AEA to inactivation, to build a pharmacokinetic model to better characterize the pharmacokinetics of peg-asparaginase and assess whether an increased clearance relates to subsequent inactivation. The study analyzed 1631 real-time AEA samples from 253 patients, identifying inactivation in 18.2% of the patients. This inactivation presented as mild allergy (28.3%), severe allergy (50.0%), or silent inactivation (21.7%). A pharmacokinetic transit compartment model was used to describe AEA-time profiles, revealing that 93% of patients with inactivation exhibited prior increased clearance, whereas 86% of patients without hypersensitivity maintained stable clearance throughout asparaginase treatment. These findings enable prediction of inactivation and options for either dose increments or a shift to alternative asparaginase formulations to optimize ALL treatment strategies.

Bridging responses to a human telomerase reverse transcriptase-based peptide cancer vaccine candidate in a mechanism-based model.

Therapeutic cancer vaccines are novel immuno-therapeutics, aiming to improve clinical outcomes with other immunotherapies. However, obstacles to their successful clinical development remain, which model-informed drug development approaches may address. UV1 is a telomerase based therapeutic cancer vaccine candidate being investigated in phase I clinical trials for multiple indications. We developed a mechanism-based model structure, using a nonlinear mixed-effects modeling techniques, based on longitudinal tumor sizes (sum of the longest diameters, SLD), UV1-specific immunological assessment (stimulation index, SI) and overall survival (OS) data obtained from a UV1 phase I trial including non-small cell lung cancer (NSCLC) patients and a phase I/IIa trial including malignant melanoma (MM) patients. The final structure comprised a mechanistic tumor growth dynamics (TGD) model, a model describing the probability of observing a UV1-specific immune response (SI ≥ 3) and a time-to-event model for OS. The mechanistic TGD model accounted for the interplay between the vaccine peptides, immune system and tumor. The model-predicted UV1-specific effector CD4+ T cells induced tumor shrinkage with half-lives of 103 and 154 days in NSCLC and MM patients, respectively. The probability of observing a UV1-specific immune response was mainly driven by the model-predicted UV1-specific effector and memory CD4+ T cells. A high baseline SLD and a high relative increase from nadir were identified as main predictors for a reduced OS in NSCLC and MM patients, respectively. Our model predictions highlighted that additional maintenance doses, i.e. UV1 administration for longer periods, may result in more sustained tumor size shrinkage.

Multistate modeling for survival analysis in critically ill patients treated with meropenem.

Appropriate antibiotic dosing to ensure early and sufficient target attainment is crucial for improving clinical outcome in critically ill patients. Parametric survival analysis is a preferred modeling method to quantify time-varying antibiotic exposure - response effects, whereas bias may be introduced in hazard functions and survival functions when competing events occur. This study investigated predictors of in-hospital mortality in critically ill patients treated with meropenem by pharmacometric multistate modeling. A multistate model comprising five states (ongoing meropenem treatment, other antibiotic treatment, antibiotic treatment termination, discharge, and death) was developed to capture the transitions in a cohort of 577 critically ill patients treated with meropenem. Various factors were investigated as potential predictors of the transitions, including patient demographics, creatinine clearance calculated by Cockcroft-Gault equation (CLCRCG ), time that unbound concentrations exceed the minimum inhibitory concentration (fT>MIC ), and microbiology-related measures. The probabilities to transit to other states from ongoing meropenem treatment increased over time. A 10 mL/min decrease in CLCRCG was found to elevate the hazard of transitioning from states of ongoing meropenem treatment and antibiotic treatment termination to the death state by 18%. The attainment of 100% fT>MIC significantly increased the transition rate from ongoing meropenem treatment to antibiotic treatment termination (by 9.7%), and was associated with improved survival outcome. The multistate model prospectively assessed predictors of death and can serve as a useful tool for survival analysis in different infection scenarios, particularly when competing risks are present.

Synergy of polymyxin B and minocycline against KPC-3- and OXA-48-producing Klebsiella pneumoniae in dynamic time-kill experiments: agreement with in silico predictions.

OBJECTIVES: Combination therapy is often used for carbapenem-resistant Gram-negative bacteria. We previously demonstrated synergy of polymyxin B and minocycline against carbapenem-resistant Klebsiella pneumoniae in static time-kill experiments and developed an in silico pharmacokinetic/pharmacodynamic (PK/PD) model. The present study assessed the synergistic potential of this antibiotic combination in dynamic experiments. METHODS: Two clinical K. pneumoniae isolates producing KPC-3 and OXA-48 (polymyxin B MICs 0.5 and 8 mg/L, and minocycline MICs 1 and 8 mg/L, respectively) were included. Activities of the single drugs and the combination were assessed in 72 h dynamic time-kill experiments mimicking patient pharmacokinetics. Population analysis was performed every 12 h using plates containing antibiotics at 4× and 8× MIC. WGS was applied to reveal resistance genes and mutations. RESULTS: The combination showed synergistic and bactericidal effects against the KPC-3-producing strain from 12 h onwards. Subpopulations with decreased susceptibility to polymyxin B were frequently detected after single-drug exposures but not with the combination. Against the OXA-48-producing strain, synergy was observed between 4 and 8 h and was followed by regrowth. Subpopulations with decreased susceptibility to polymyxin B and minocycline were detected throughout experiments. For both strains, the observed antibacterial activities showed overall agreement with the in silico predictions. CONCLUSIONS: Polymyxin B and minocycline in combination showed synergistic effects, mainly against the KPC-3-producing K. pneumoniae. The agreement between the experimental results and in silico predictions supports the use of PK/PD models based on static time-kill data to predict the activity of antibiotic combinations at dynamic drug concentrations.

Model-Informed Precision Dosing to Reduce Vincristine-Induced Peripheral Neuropathy in Pediatric Patients: A Pharmacokinetic and Pharmacodynamic Modeling and Simulation Analysis.

BACKGROUND: Vincristine-induced peripheral neuropathy (VIPN) is a common adverse effect of vincristine, a drug often used in pediatric oncology. Previous studies demonstrated large inter- and intrapatient variability in vincristine pharmacokinetics (PK). Model-informed precision dosing (MIPD) can be applied to calculate patient exposure and individualize dosing using therapeutic drug monitoring (TDM) measurements. This study set out to investigate the PK/pharmacodynamic (PKPD) relationship of VIPN and determine the utility of MIPD to support clinical decisions regarding dose selection and individualization. METHODS: Data from 35 pediatric patients were utilized to quantify the relationship between vincristine dose, exposure and the development of VIPN. Measurements of vincristine exposure and VIPN (Common Terminology Criteria for Adverse Events [CTCAE]) were available at baseline and for each subsequent dosing occasions (1-5). A PK and PKPD analysis was performed to assess the inter- and intraindividual variability in vincristine exposure and VIPN over time. In silico trials were performed to portray the utility of vincristine MIPD in pediatric subpopulations with a certain age, weight and cytochrome P450 (CYP) 3A5 genotype distribution. RESULTS: A two-compartmental model with linear PK provided a good description of the vincristine exposure data. Clearance and distribution parameters were related to bodyweight through allometric scaling. A proportional odds model with Markovian elements described the incidence of Grades 0, 1 and ≥ 2 VIPN overdosing occasions. Vincristine area under the curve (AUC) was the most significant exposure metric related to the development of VIPN, where an AUC of 50 ng⋅h/mL was estimated to be related to an average VIPN probability of 40% over five dosing occasions. The incidence of Grade ≥ 2 VIPN reduced from 62.1 to 53.9% for MIPD-based dosing compared with body surface area (BSA)-based dosing in patients. Dose decreases occurred in 81.4% of patients with MIPD (vs. 86.4% for standard dosing) and dose increments were performed in 33.4% of patients (no dose increments allowed for standard dosing). CONCLUSIONS: The PK and PKPD analysis supports the use of MIPD to guide clinical dose decisions and reduce the incidence of VIPN. The current work can be used to support decisions with respect to dose selection and dose individualization in children receiving vincristine.

A model-based approach leveraging in vitro data to support dose selection from the outset: A framework for bispecific antibodies in immuno-oncology.

FAP-4-1BBL is a bispecific antibody exerting 4-1BB-associated T-cell activation only while simultaneously bound to the fibroblast activation protein (FAP) receptor, expressed on the surface of cancer-associated fibroblasts. The trimeric complex formed when FAP-4-1BBL is simultaneously bound to FAP and 4-1BB represents a promising mechanism to achieve tumor-specific 4-1BB stimulation. We integrated in vitro data with mathematical modeling to characterize the pharmacology of FAP-4-1BBL as a function of trimeric complex formation when combined with the T-cell engager cibisatamab. This relationship was used to prospectively predict a range of clinical doses where trimeric complex formation is expected to be at its maximum. Depending on the dosing schedule and FAP-4-1BBL plasma: tumor distribution, doses between 2 and 145 mg could lead to maximum trimeric complex formation in the clinic. Due to the expected variability in both pharmacokinetic and FAP expression in the patient population, we predict that detecting a clear dose-response relationship would remain difficult without a large number of patients per dose level, highlighting that mathematical modeling techniques based on in vitro data could aid dose selection.

Model-based assessment of neutrophil-mediated phagocytosis and digestion of bacteria across in vitro and in vivo studies.

Neutrophil granulocytes are key components of the host response against pathogens, and severe neutropenia, with neutrophil counts below 0.5 × 106 cells/mL, renders patients increasingly vulnerable to infections. Published in vitro (n = 7) and in vivo (n = 5) studies with time-course information on bacterial and neutrophil counts were digitized to characterize the kinetics of neutrophil-mediated bacterial killing and inform on the immune systems' contribution to the clearance of bacterial infections. A mathematical model for the in vitro dynamics of bacteria and the kinetics of neutrophil-mediated phagocytosis and digestion was developed, which was extended to in vivo studies in immune-competent and immune-compromised mice. Neutrophil-mediated bacterial killing was described by two first-order processes-phagocytosis and digestion-scaled by neutrophil concentration, where 50% of the maximum was achieved at neutrophil counts of 1.19 × 106 cells/mL (phagocytosis) and 6.55 × 106 cells/mL (digestion). The process efficiencies diminished as the phagocytosed bacteria to total neutrophils ratio increased (with 50% reduction at a ratio of 3.41). Neutrophil in vivo dynamics were captured through the characterization of myelosuppressive drug effects and postinoculation neutrophil influx into lungs and by system differences (27% bacterial growth and 9.3% maximum capacity, compared with in vitro estimates). Predictions showed how the therapeutically induced reduction of neutrophil counts enabled bacterial growth, especially when falling below 0.5 × 106 cells/mL, whereas control individuals could deal with all tested bacterial burdens (up to 109 colony forming units/g lung). The model-based characterization of neutrophil-mediated bacterial killing simultaneously predicted data across in vitro and in vivo studies and may be used to inform the capacity of host-response at the individual level.

Predicting cytokine kinetics during sepsis; a modelling framework from a porcine sepsis model with live Escherichia coli.

BACKGROUND: Describing the kinetics of cytokines involved as biomarkers of sepsis progression could help to optimise interventions in septic patients. This work aimed to quantitively characterise the cytokine kinetics upon exposure to live E. coli by developing an in silico model, and to explore predicted cytokine kinetics at different bacterial exposure scenarios. METHODS: Data from published in vivo studies using a porcine sepsis model were analysed. A model describing the time courses of bacterial dynamics, endotoxin (ETX) release, and the kinetics of TNF and IL-6 was developed. The model structure was extended from a published model that quantifies the ETX-cytokines relationship. An external model evaluation was conducted by applying the model to literature data. Model simulations were performed to explore the sensitivity of the host response towards differences in the input rate of bacteria, while keeping the total bacterial burden constant. RESULTS: The analysis included 645 observations from 30 animals. The blood bacterial count was well described by a one-compartment model with linear elimination. A scaling factor was estimated to quantify the ETX release by bacteria. The model successfully described the profiles of TNF, and IL-6 without a need to modify the ETX-cytokines model structure. The kinetics of TNF, and IL-6 in the external datasets were well predicted. According to the simulations, the ETX tolerance development results in that low initial input rates of bacteria trigger the lowest cytokine release. CONCLUSION: The model quantitively described and predicted the cytokine kinetics triggered by E. coli exposure. The host response was found to be sensitive to the bacterial exposure rate given the same total bacterial burden.

Assessment of ibrutinib scheduling on leukocyte, lymph node size and blood pressure dynamics in chronic lymphocytic leukemia through pharmacokinetic-pharmacodynamic models.

Ibrutinib is a Bruton tyrosine kinase (Btk) inhibitor for treating chronic lymphocytic leukemia (CLL). It has also been associated with hypertension. The optimal dosing schedule for mitigating this adverse effect is currently under discussion. A quantification of relationships between systemic ibrutinib exposure and efficacy (i.e., leukocyte count and sum of the product of perpendicular diameters [SPD] of lymph nodes) and hypertension toxicity (i.e., blood pressure), and their association with overall survival is needed. Here, we present a semi-mechanistic pharmacokinetic-pharmacodynamic modeling framework to characterize such relationships and facilitate dose optimization. Data from a phase Ib/II study were used, including ibrutinib plasma concentrations to derive daily 0-24-h area under the concentration-time curve, leukocyte count, SPD, survival, and blood pressure measurements. A nonlinear mixed effects modeling approach was applied, considering ibrutinib's pharmacological action and CLL cell dynamics. The final framework included (i) an integrated model for SPD and leukocytes consisting of four CLL cell subpopulations with ibrutinib inhibiting phosphorylated Btk production, (ii) a turnover model in which ibrutinib stimulates an increase in blood pressure, and (iii) a competing risk model for dropout and death. Simulations predicted that the approved dosing schedule had a slightly higher efficacy (24-month, progression-free survival [PFS] 98%) than de-escalation schedules (24-month, average PFS ≈ 97%); the latter had, on average, ≈20% lower proportions of patients with hypertension. The developed modeling framework offers an improved understanding of the relationships among ibrutinib exposure, efficacy and toxicity biomarkers. This framework can serve as a platform to assess dosing schedules in a biologically plausible manner.

A multistate modeling and simulation framework to learn dose-response of oncology drugs: Application to bintrafusp alfa in non-small cell lung cancer.

The dose/exposure-efficacy analyses are often conducted separately for oncology end points like best overall response, progression-free survival (PFS) and overall survival (OS). Multistate models offer to bridge these dose-end point relationships by describing transitions and transition times from enrollment to response, progression, and death, and evaluating transition-specific dose effects. This study aims to apply the multistate pharmacometric modeling and simulation framework in a dose optimization setting of bintrafusp alfa, a fusion protein targeting TGF-ß and PD-L1. A multistate model with six states (stable disease [SD], response, progression, unknown, dropout, and death) was developed to describe the totality of endpoints data (time to response, PFS, and OS) of 80 patients with non-small cell lung cancer receiving 500 or 1200 mg of bintrafusp alfa. Besides dose, evaluated predictor of transitions include time, demographics, premedication, disease factors, individual clearance derived from a pharmacokinetic model, and tumor dynamic metrics observed or derived from tumor size model. We found that probabilities of progression and death upon progression decreased over time since enrollment. Patients with metastasis at baseline had a higher probability to progress than patients without metastasis had. Despite dose failed to be statistically significant for any individual transition, the combined effect quantified through a model with dose-specific transition estimates was still informative. Simulations predicted a 69.2% probability of at least 1 month longer, and, 55.6% probability of at least 2-months longer median OS from the 1200 mg compared to the 500 mg dose, supporting the selection of 1200 mg for future studies.

Prevalence and Clinical Consequences of Colistin Heteroresistance and Evolution into Full Resistance in Carbapenem-Resistant Acinetobacter baumannii.

Colistin heteroresistance (HR) refers to a bacterial population comprised of several subpopulations with different levels of resistance to colistin. In this study, we discuss the classic form of HR, in which a resistant subpopulation exists within a predominantly susceptible population. We investigated the prevalence of colistin HR and its evolution into full resistance among 173 clinical carbapenem-resistant Acinetobacter baumannii isolates and examined the effect of HR on clinical outcomes. To determine HR, we performed population analysis profiling. Our results showed a high prevalence of HR (67.1%). To examine evolution of HR strains into full resistance, the HR strains were grown in colistin-containing broth, transferred onto colistin-containing plates, and colonies on these plates were transferred into colistin-free broth. Many of the HR strains (80.2%) evolved into full resistance, 17.2% reverted to HR, and 2.6% were borderline. We used logistic regression to compare 14-day clinical failure and 14-day mortality between patients infected by HR versus susceptible non-HR carbapenem-resistant A. baumannii. In the subgroup of patients with bacteremia, HR was significantly associated with 14-day mortality. IMPORTANCE To our knowledge, this is the first large-scale study to report on HR in Gram-negative bacteria. We described the prevalence of colistin HR in a large sample of carbapenem-resistant A. baumannii isolates, the evolution of many colistin HR isolates to a resistant phenotype following colistin exposure and withdrawal, and the clinical consequences of colistin HR. We found a high prevalence of HR among clinical carbapenem-resistant A. baumannii isolates; most evolved into a resistant phenotype following colistin exposure and withdrawal. In patients treated with colistin, evolution of HR A. baumannii into full resistance could lead to higher rates of treatment failure and contribute to the reservoir of colistin-resistant pathogens in health care settings.

Multistate Pharmacometric Model to Define the Impact of Second-Line Immunotherapies on the Survival Outcome of the IMpower131 Study.

Overall survival is defined as the time since randomization into the clinical trial to event of death or censor (end of trial or follow-up), and is considered to be the most reliable cancer end point. However, the introduction of second-line treatment after disease progression could influence survival and be considered a confounding factor. The aim of the current study was to set up a multistate model framework, using data from the IMpower131 study, to investigate the influence of second-line immunotherapies on overall survival analysis. The model adequately described the transitions between different states in patients with advanced squamous non-small cell lung cancer treated with or without atezolizumab plus nab-paclitaxel and carboplatin, and characterized the survival data. High PD-L1 expression at baseline was associated with a decreased hazard of progression, while the presence of liver metastasis at baseline was indicative of a high risk of disease progression after initial response. The hazard of death after progression was lower for participants who had longer treatment response, i.e., longer time to progression. The simulations based on the final multistate model showed that the addition of atezolizumab to the nab-paclitaxel and carboplatin regimen had significant improvement in the patients' survival (hazard ratio = 0.75, 95% prediction interval: 0.61-0.90 favoring the atezolizumab + nab-paclitaxel and carboplatin arm). The developed modeling approach can be applied to other cancer types and therapies to provide a better understanding of efficacy of drug and characterizing different states, and investigate the benefit of primary therapy in survival while accounting for the switch to alternative treatment in the case of disease progression.

Optimization of blood pressure measurement practices for pharmacodynamic analyses of tyrosine-kinase inhibitors.

Blood pressure measurements form a critical component of adverse event monitoring for tyrosine kinase inhibitors, but might also serve as a biomarker for dose titrations. This study explored the impact of various sources of within-individual variation on blood pressure readings to improve measurement practices and evaluated the utility for individual- and population-level dose selection. A pharmacokinetic-pharmacodynamic modeling framework was created to describe circadian blood pressure changes, inter- and intra-day variability, changes from dipper to non-dipper profiles, and the relationship between drug exposure and blood pressure changes over time. The framework was used to quantitatively evaluate the influence of physiological and pharmacological aspects on blood pressure measurements, as well as to compare measurement techniques, including office-based, home-based, and ambulatory 24-h blood pressure readings. Circadian changes, as well as random intra-day and inter-day variability, were found to be the largest sources of within-individual variation in blood pressure. Office-based and ambulatory 24-h measurements gave rise to potential bias (>5 mmHg), which was mitigated by model-based estimations. Our findings suggest that 5-8 consecutive, home-based, measurements taken at a consistent time around noon, or alternatively within a limited time frame (e.g., 8.00 a.m. to 12.00 p.m. or 12.00 p.m. to 5.00 p.m.), will give rise to the most consistent blood pressure estimates. Blood pressure measurements likely do not represent a sufficiently accurate method for individual-level dose selection, but may be valuable for population-level dose identification. A user-friendly tool has been made available to allow for interactive blood pressure simulations and estimations for the investigated scenarios.

Rituximab pharmacokinetic and pharmacokinetic-pharmacodynamic evaluation based on a study in diffuse large B-cell lymphoma: Influence of tumor size on pharmacokinetic and assessment of pharmacokinetic similarity.

Dr. Reddy's Laboratories rituximab (DRL_RI; Dr. Reddy's Laboratories SA, Basel, Switzerland) is under development as a rituximab biosimilar. Study RI-01-002 (Clinical Trials Registry - India/2012/11/003129), comparing DRL_RI to the reference medicinal product (RMP) MabThera® (Roche, Grenzach-Wyhlen, Germany), demonstrated pharmacokinetic (PK) equivalence and showed comparable pharmacodynamic, efficacy, safety, and immunogenicity profiles. We used data from the same study to perform population PK and PK-pharmacodynamic analyses: first exploring possible factors influencing the PK similarity assessment between products and then performing simulations to investigate the impact of tumor size on rituximab PK. Nonlinear mixed-effects models for PK, tumor size, tumor size-PK, and tumor response were developed independently. The final PK model included drug product as a dose-scaling parameter and predicted a 6.75% higher dose reaching the system in RMP-treated patients. However, when tumor size was included in the tumor size-PK model, the drug product effect was no longer observed. The model rather indicated that patients with larger tumor size have higher clearance. Further simulations confirmed that higher baseline tumor size is associated to slightly lower rituximab exposure. Tumor response, described by a continuous-time Markov model, did not differ between drug products. Both had higher effects during the first 20 weeks of treatment. Also, the model described a subpopulation of nonresponders to treatment (42%) with faster transitions to a worse state. The different rituximab exposure initially detected between drug products (6.75%) was shown using PK/PK-pharmacodynamic analysis to be attributed to a tumor size imbalance between treatment groups. PK/PK-pharmacodynamic analyses may contribute to PK similarity assessments.

Bayesian forecasting of tumor size metrics and overall survival.

The tumor size ratio (TSR), time-to-tumor growth (TTG), and tumor growth rate (kG) are frequently suggested as model-based predictors of overall survival (OS) for different types of tumors. When the tumor metrics are applied in forecasting of the outcome for individual patients at an early stage, the tumor data might be sparse resulting in imprecise prediction. This simulation study aimed to investigate how the tumor follow-up data and estimation approaches influence the accuracy in the tumor size metrics and the predicted hazard of death for individual patients. Longitudinal tumor size and OS data were simulated using tumor growth inhibition and Weibull distribution models, respectively. Based on the model and increasing measurement durations, the accuracy (defined as 80-125% of the simulated "true" value) in individual metrics and hazard was computed. TSR week 6 (TSRw6) accuracy was adequate for 91% of the patients when tumor size was measured up to 12 weeks. For TTG and kG metrics, the highest accuracy observed was lower (43 and 77%, respectively) and occurred later (42 and 60 weeks, respectively). The simultaneous (joint) and sequential estimation approaches resulted in similar accuracies, however, in general, the sequential approach where individual tumor size parameters are fixed, demonstrated inferior estimation properties. The TSRw6 and the model-predicted tumor time course (absolute or relative change) had better forecasting properties than TTG or kG. The population pharmacokinetic (PK) parameters and data approach performed similarly or better than the simultaneous approach and had a better accuracy in estimating individuals' hazard of death than the individual PK parameters method.

Preclinical InVivo Data Integrated in a Modeling Network Informs a Refined Clinical Strategy for a CD3 T-Cell Bispecific in Combination with Anti-PD-L1.

TYRP1-TCB is a CD3 T-cell bispecific (CD3-TCB) antibody for the treatment of advanced melanoma. A tumor growth inhibition (TGI) model was developed using mouse xenograft data with TYRP1-TCB monotherapy or TYRP1-TCB plus anti-PD-L1 combination. The model was translated to humans to inform a refined clinical strategy. From xenograft mouse data, we estimated an EC50 of 0.345 mg/L for TYRP1-TCB, close to what was observed in vitro using the same tumor cell line. The model showed that, though increasing the dose of TYRP1-TCB in monotherapy delays the time to tumor regrowth and promotes higher tumor cell killing, it also induces a faster rate of tumor regrowth. Combination with anti-PD-L1 extended the time to tumor regrowth by 25% while also decreasing the tumor regrowth rate by 69% compared to the same dose of TYRP1-TCB alone. The model translation to humans predicts that if patients' tumors were scanned every 6 weeks, only 46% of the monotherapy responders would be detected even at a TYRP1-TCB dose resulting in exposures above the EC90. However, combination of TYRP1-TCB and anti-PD-L1 in the clinic is predicted to more than double the overall response rate (ORR), duration of response (DoR) and progression-free survival (PFS) compared to TYRP1-TCB monotherapy. As a result, it is highly recommended to consider development of CD3-TCBs as part of a combination therapy from the outset, without the need to escalate the CD3-TCB up to the Maximum Tolerated Dose (MTD) in monotherapy and without gating the combination only on RECIST-derived efficacy metrics.

Expert workshop summary: Advancing toward a standardized murine model to evaluate treatments for antimicrobial resistance lung infections.

The rise in antimicrobial resistance (AMR), and increase in treatment-refractory AMR infections, generates an urgent need to accelerate the discovery and development of novel anti-infectives. Preclinical animal models play a crucial role in assessing the efficacy of novel drugs, informing human dosing regimens and progressing drug candidates into the clinic. The Innovative Medicines Initiative-funded "Collaboration for prevention and treatment of MDR bacterial infections" (COMBINE) consortium is establishing a validated and globally harmonized preclinical model to increase reproducibility and more reliably translate results from animals to humans. Toward this goal, in April 2021, COMBINE organized the expert workshop "Advancing toward a standardized murine model to evaluate treatments for AMR lung infections". This workshop explored the conduct and interpretation of mouse infection models, with presentations on PK/PD and efficacy studies of small molecule antibiotics, combination treatments (ß-lactam/ß-lactamase inhibitor), bacteriophage therapy, monoclonal antibodies and iron sequestering molecules, with a focus on the major Gram-negative AMR respiratory pathogens Pseudomonas aeruginosa, Klebsiella pneumoniae and Acinetobacter baumannii. Here we summarize the factors of variability that we identified in murine lung infection models used for antimicrobial efficacy testing, as well as the workshop presentations, panel discussions and the survey results for the harmonization of key experimental parameters. The resulting recommendations for standard design parameters are presented in this document and will provide the basis for the development of a harmonized and bench-marked efficacy studies in preclinical murine pneumonia model.